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Abstract 

We present a comprehensive analysis of the recent data from the E866 ex- 
periment at Fermilab on Drell-Yan production in pD and pp collisions, which 
indicates a non-trivial x-dependence for the asymmetry between u and d quark 
distributions in the proton. The relatively fast decrease of the asymmetry at 
large x suggests the important role played by the chiral structure of the nu- 
cleon, in particular the irN and irA components of the nucleon wave function. 
At small x the data require an additional non-chiral component, which may 
be attributed to the Pauli exclusion principle, as first suggested by Field and 
Feynman. 

I. INTRODUCTION 

The recent Drell-Yan experiment by the E866/NuSea Collaboration at Fermilab [|lj pro- 
vides the best information yet on the detailed structure of the light antiquark sea of the 
proton. Although previous experiments by the New Muon Collaboration || on the differ- 
ence between the proton and neutron structure functions established that an asymmetry in 
the sea exists, they yielded direct information only on the first moment of the antiquark 
asymmetry. The earlier Drell-Yan experiment by the NA51 Collaboration at CERN || 
measured the up and down antiquark ratio, though at zero rapidity |4j], but in order to 
improve the statistical accuracy the data were binned to a single value of Bjorken-x. The 
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E866 experiment, on the other hand, has for the first time mapped out the shape of the d/u 
ratio over a large range of x, 0.02 < x < 0.345. 

The relatively large asymmetry found in these experiments implies the presence of non- 
trivial dynamics in the proton sea which does not have a perturbative QCD origin. From the 
symmetry properties of QCD, we know that one source of non-perturbative quark- ant iquark 
pairs is the pion cloud associated with dynamical chiral symmetry breaking. The novel 
and unexpected feature of the E866 data is that the d/u asymmetry peaks at rather small 
values of x, and drops quite rapidly with increasing x, approximately like (1 — x) n with 
n ~ 10. This behavior had not been anticipated in global parameterizations of data, and 
is softer than pion cloud models of the nucleon would generally predict. While the former 
is just an artifact of an extrapolation of a parametric function into an unmeasured region, 
with no physics implications, the reason for the latter is the rather hard valence antiquark 
distribution in the pion, q£ ~ (1 — x) n with n ~ 1. Although this x dependence is softened 
somewhat by the probability of finding the pion in the nucleon in the first place, the resulting 
convoluted distribution still does not die off as rapidly, at large x, as the new E866 data 
would suggest. 

Our analysis suggests that a quantitative description of the entire region of x covered in 
the experiment requires a delicate balance between several competing mechanisms, which 
leads us to speculate about the possibility of a two-phase picture of the non-perturbative 
sea of the nucleon. At larger x, the dynamics of the pion cloud of the nucleon come to the 
fore, with deep-inelastic scattering from the nN component of the nucleon wave function 
providing the bulk of the d — u asymmetry. On the other hand, the important role of the 
A isobar in nuclear physics has been known for a long time, and it proves to be of some 
importance in this case too. If a part of the 7rA distribution happens to be harder than 
the irN, there would be cancellation of some of the d excess at large x. Such a piece does 
indeed arise in the light-cone formulation of the meson-cloud model, provided the ttNA form 
factor, parameterized for example by a dipole form, is harder than the nNN form factor 
H — something which is consistent with the measured difference between the nucleon and 
NA transition axial form factors. 

To be consistent with the measured sum, u + d, it is known that both the 7iNN and nNA 
form factors need to be relatively soft 0. The best fit, within the pion cloud framework, to 
the data on both the sum and difference of u and d at large x accounts for around half of the 
integrated asymmetry, leaving room for possible other, non-pionic, mechanisms to provide 
the missing strength at smaller x. (While the usual discussions of the pionic contribution 
focus on the valence quarks in the pion, there is also some theoretical argument for a d — u 
asymmetry in the sea of the pion. We also estimate how this might affect the analysis.) 

It should be noted that, aside from the flavor- asymmetric nN and 7rA components of the 
nucleon, there is no a priori reason why the 'bare' (non-pion-dressed) nucleon state itself 
cannot have an intrinsic asymmetric sea associated with it. In fact, this is actually what is 
expected from the Pauli exclusion principle, as anticipated long ago by Field and Feynman 
pj on the simple basis that the u and d valence quark sectors are unequally populated 
in the proton ground state. Although more difficult to estimate model-independently, the 
contribution to the d — u difference from antisymmetrization has been calculated within a 
non-perturbative model of the nucleon ||. Along the lines of the model estimates, we find 
that the effects of antisymmetrization are most relevant at small x, with normalization such 
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that they can account for most of the remaining half of the integrated asymmetry. Indeed, 
our analysis suggests that the best fit to the E866 data is obtained when both mechanisms 
play a role, consistent with the conclusions of earlier analyses PJTOf of the NMC data on the 
proton-neutron structure function difference. 

In the following, we firstly outline in Section II the experimental status of the d/u 
asymmetry, including a comparison of the Drell-Yan data with deep-inelastic scattering 
data. In Section III the question of the possible origin of the asymmetry is addressed, in 
the form of the chiral structure of the nucleon and the associated pion cloud, while Section 
IV deals with the role of the Pauli exclusion principle as a source of flavor asymmetry. In 
Section V we summarize our conclusions. 



II. THE LIGHT ANTIQUARK ASYMMETRY 

The E866 /NuSea Collaboration measured /i + /U~ Drell-Yan pairs produced in pp and pD 
collisions. In the parton model the Drell-Yan cross section is proportional to: 

o** oc £ e\ (f{ Xl ) q h (x 2 ) + <f{x x ) q\x 2 )) , (1) 
<? 

where h = p or D, and x\ and x 2 are the light-cone momentum fractions carried by partons 
in the projectile and target hadron, respectively. 

Assuming that the deuteron is composed of two bound nucleons, and utilizing isospin 
symmetry {u p = d n , etc.), in the limit X\ ^> x 2 (in which q{x\) <C q(xi)) the ratio of the 
deuteron to proton cross sections can be written: 
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4 + d(x 1 )/u(xi) ^ 



4 + d{x\)/u{x\) ■ d(x 2 )/u(x 2 ) 



where q is the antiquark distribution in the bound proton. Neglecting relativistic and nucleon 
off-shell effects, this can be approximated by a convolution of the antiquark distribution in 
the proton with the proton distribution function in the deuteron |TT],IT2|, 



~ r dz 

q(x) « / —f N/D (z) q(x/z), (3) 

J x Z 

where Jn/dI 2 ) is the distribution of nucleons in the deuteron with light-cone momentum 



fraction z. In practice we use the function /jv/d from Ref. where it is given in terms of 
a realistic deuteron wave function that has been constrained to reproduce the static deuteron 
properties and nucleon-nucleon phase shifts. 

In the absence of nuclear effects, q = q, one would have: 

<J pD = 1 L | d{x 2 ) \ A + d{x l )/u{x l ) 

Icjpp 2\ u(x 2 )J A + d(x 1 )/u(x 1 ) ■ d(x 2 )/u(x 2 ) ' 

so that the ratio would be unity if d — u. On the other hand, we know that nuclear 



shadowing exists in the deuteron at small x (see [0 and references within), and at large 



x nuclear binding and Fermi motion effects come into prominence W%- Since the bulk of 
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the effect is observed outside the very small x region, shadowing will not be important 
in these data. However, Fermi smearing is potentially more relevant, as what enters in 
the parentheses in Eq.(|2|) is the ratio of smeared to unsmeared sea quark distributions, for 
which smearing effects should come into play at much smaller x than for the total structure 



function [IT] . In Fig. 1 we show the ratio of the antiquark distributions in a proton bound 
in the deuteron to that in a free proton, q/q. Parameterizing the antiquarks for illustration 
purposes at large a; by a simple q ~ (1 — x) n form, with n = 5, 7 and 10, the ratio is seen to 
rise rapidly above x ~ 0.4, though in the measured region it only deviates from unity by a 
few percent. From this one can conclude that Eq.(^) should be a reasonable approximation 
to Eq.(@). 

In the extreme large-^i limit, where d(x\) <C u(xi), the cross section ratio would directly 
give d/u: 
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(zi)<u(a;i) 



For the E866 data the criterion x\ ^> x<i is not always satisfied, however, so that Eq.@ gives 
only an indication of the sensitivity of the Drell-Yan cross section to d/u, and in practice 
the full parton model cross section is used together with an iterative procedure in which the 
valence and total sea distributions are assumed known and the extracted d/u ratio adjusted 
to fit the data Q. The resulting d/u ratio is shown in Fig. 2, where the asymmetry is found 
to peak at relatively small x, x ~ 0.15, dropping rapidly to unity by x ~ 0.3. Also shown for 
comparison is the NA51 data point [0], extracted from the Drell-Yan pp and pn asymmetry 
at x = 0.18, which lies slightly above the E866 data. The two fits are from the CTEQ4 [15| 



parameterization (dashed line), which has the asymmetry extending out to larger x, and the 
more recent MRS98 analysis |16| (dotted line), which included the E866 data. 

From the d/u ratio the E866/NuSea Collaboration further extract the difference, A 
(I a . assuming the sum d + u from the CTEQ4 fit, 

\d/ U +l J V / CTEQ4 V ; 

The resulting points are shown in Fig. 3, again in comparison with the CTEQ4 and MRS98 
parameterizations. It should be noted, however, that the sea quark distribution in these 
fits is not very well determined at large x. In particular, a larger total sea at x ~ 0.2-0.3 
would result in a larger asymmetry A in the region just where the E866 data appear to drop 
rapidly to zero. Further data from the E866 experiment on the total antiquark distribution 
at large x should help to clarify the issue. The Drell-Yan data can also be compared with 
the proton-neutronF] structure function difference measured previously by the New Muon 
Collaboration [0], if one assumes that the valence quark distributions in the proton are 
known. In this case the antiquark asymmetry can be written: 



1 Note that the NMC extracted the neutron structure function from proton and deuteron data 
assuming F£ = F% — F%. 
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,h.r) - U{X) = \ (U V (X) - dy(x)) - |- (Ff (*) - F 2 ») NMC . (7) 



The d — u difference extracted from the NMC data is shown in Fig. 3 using both the CTEQ4 
(open circles) and MRS98 (diamonds) fits to the valence quark distributions. The NMC 
values appear to lie consistently above the E866 data for x above ~ 0.1, although there 
is some sensitivity to the choice of valence quark parameterization. Not surprisingly, the 
E866 integrated value (after extrapolating down to x = and up to x = 1) is found to be: 
A E866 = j^dx(d-u) = 0.100 ± 0.018 13, somewhat smaller than the NMC value, which 



is A NMC = 0.148 ± 0.039 0, although still consistent within errors. This difference will be 
further enhanced if one corrects the NMC data for shadowing in the deuteron, omission of 



which underestimates the violation of the Gottfried sum rule |L3 . 

The observation of a large asymmetry between u and d, now both at CERN and Fermi- 
lab, provides theorists with a challenge to better understand the internal, non-perturbative 
structure of the nucleon, as the asymmetry due to perturbative effects is known to be very 
small |18|]. In the next section we consider the possible origin of this asymmetry, in the form 



of the non-perturbative chiral structure of the nucleon. 



III. CHIRAL SYMMETRY AND THE MESON CLOUD 



The simplest and most obvious source of a non-perturbative asymmetry in the light 
quark sea is the chiral structure of QCD. From numerous studies in low energy physics, 
including chiral perturbation theory, pions are known to play a crucial role in the structure 
and dynamics of the nucleon. However, there is no reason why the long-range tail of the 
nucleon should not also play a role at higher energies. This was first alluded to by Sullivan 
T9| , who argued that deep- inelastic scattering from the pion cloud of the nucleon is a scaling 



contribution to the nucleon structure function. Indeed, expectations for the ratio of nuclear 
to nucleon structure functions, based on arguments that nuclear scales are much smaller 
than typical deep- inelastic scales, and therefore irrelevant, were proved to be dramatically 



wrong by the observation of the nuclear EMC effect p0 

As pointed out by Thomas |J, if the proton's wave function contains an explicit 7i + n 
Fock state component, a deep-inelastic probe scattering from the virtual 7r + , which contains 
a valence d quark, will automatically lead to ad excess in the proton. This is the essential 
physical idea behind these expectations, and has been used to address not only the d/u 
asymmetry [p|,p|, p!0| , pT| -|27[] , but also SU(3) flavor symmetry breaking in the proton sea ||, 
as well as asymmetries in the strange 



29,30]. In recent years 



p8| and heavier flavor sectors 
this picture has been refined and elaborated with inclusion of additional meson and baryon 
states 1, and constraints put on many of the model parameters by comparisons of the 



predictions of the model with other processes [31.32 



The basic hypothesis of the meson cloud model is that the physical nucleon state can be 
expanded (in an infinite momentum frame (IMF) and in the one-meson approximation) in 
a series involving bare nucleon and two-particle, meson-baryon states. The essential ingre- 
dients are the meson-baryon distribution functions, /mb(|/), which give the probability to 
find a meson, M, in the nucleon carrying a fraction y of the nucleon's light-cone momentum. 
As discussed at length in the literature, for these functions to have the correct probabilis- 
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tic interpretation in the IMF, they must be related to the distributions of baryons in the 
nucleon, f B Nf(y), via: 



fMB(y) = f. 



BM\ 



y)- 



(8) 



This constraint can be verified easily in the IMF, but is not entirely clear in covariant 
formulations (see Refs. pj23,27| for further discussion on this point). The IMF treatment 
has the additional advantage that the meson and baryon are on-mass-shell and so one has no 
ambiguities associated with the possible off-mass- shell behavior of their structure functions 
that are encountered in the covariant treatments. 

The contribution to the antiquark distribution in the proton, 5^ MB ^q, can then be written 
in the IMF as a convolution of the meson distribution function and the antiquark distribution 
in the (on-mass-shell) pion: 



Jx y 



(y)q M (x/y). 



(9) 



Note that this is the leading contribution to the antiquark distribution, and is independent 
of the model of the bare nucleon states. 

In earlier studies it was found, not surprisingly, that pions are the most important mesons, 
and that the dominant contributions are those associated with the ttN component of the 
proton's wave function. The distribution of pions with a recoiling nucleon is given by 

ra^ 



3g 



Inn 



16tt 2 



^tv(^tv) (k 2 T + y 2 M 2 



o (l — y)y (M 2 



■>ttN 



i-y 



(10) 



so that f^+n = 2f n o p = (2/3) f n N for the respective charge states. The invariant mass squared 
of the iiN system is given by s n N = (kj< + m 2 T )/y+ (k 2 ^ + M 2 )/(l — y), and for the functional 
form of the ttNN vertex form factor J-^n{ s -kn) we take a simple dipole parameterization: 



F kn{.S-kN) 



Kn + M 2 



^ttN + S tN , 



(11) 



normalized so that the coupling constant g n NN has its standard value (= 13.07) at the pole 
(^(M 2 ) = 1). The symmetry relations (§) are automatically guaranteed with this type of 
form factor, whereas in the earlier covariant formulations [|9t|T0,pT|-23l, with t-dependent form 
factors, this could not be achieved. Note that the E866 group also utilized the formulation 
in terms of t-dependent form factors in their recent theoretical analysis |17j of the Drell-Yan 
data. 

The antiquark distribution in the pion has been measured in ttN Drell-Yan experiments 
by the E615 Collaboration at Fermilab and by the NA10 [34] and NA3 Collabora- 
tions at CERN. These have been parameterized in next-to-leading order analyses in Refs. 
3~6| , |3~7|| . Unless stated otherwise, we use the valence part of the pion's antiquark distribution 



from Ref. [36[] throughout this analysis. 

Because the meson cloud model is a model of part of the (non-perturbative) sea, it can 
only be reliably applied to describing the non-singlet d — u distribution. In Fig. 4 we show 
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the calculated difference arising from the irN component of the proton's wave function for 
two different cut-off masses, A nN = 1 GeV (dashed) and 1.5 GeV (solid), giving average 
multiplicities (n)^ = Jq dyf v ^{y) = 13% and 26%, respectively. With the latter one has 
excellent agreement at intermediate x, x < 0.2, while somewhat overestimating the data at 
the larger x values. The excess at large x is less severe for the smaller cut-off. However, the 
strength of the contribution in that case is too small at lower x. 

To reconstruct the ratio from the calculated difference, we assume, following E866, that 
the total d + u is given by the CTEQ4 parameterization |0j|, and invert Eq.(|]). The 



resulting ratio is plotted in Fig. 5. At small x the agreement with the data for the larger 
cut-off is clearly excellent, but at larger x the calculation does not follow the downward 
trend suggested by the data, similar to the CTEQ4 parameterization in Fig. 2. 

Both Figs. 4 and 5 suggest that the excess of dover u is too strong for x ^ 0.2, and that a 
mechanism which suppresses or cancels this excess could be responsible for the behavior seen 
in the data. One way to obtain such a suppression would be if the pion structure function 
were softer. Actually, perturbative QCD suggests that the leading twist part of F£ should 
behave like (1 — x) 2 |38| at large x (see also Ref. |39[]), although the Drell-Yan data [|3"3}-|3"5 



show that it is closer to (1 — x), even with the inclusion of higher twist contributions. 

Within the pion cloud model another way that some of the cancellation can be understood 
is through the A isobar. Although the rrA component should be smaller in magnitude than 
the 7rN, if some part of the irA distribution were to be harder than the nN it would allow for 
some cancellation of the large-x excess, while preserving more of the asymmetry at smaller 
x. In fact, qualitatively such behavior is exactly what is seen in the model. 

In the IMF, the pion distribution function with a recoil A is given by P, |2~T| , |3~Tf : 

2&na r dk\ ^ A (s wA ) 



16tt 2 Jo (l-y)y (s„ A - M 2 ) 2 

^ [k 2 T + (Ma - (1 - y)M) 2 } [k 2 + (Ma + (1 - y)Mf 
6M|(l-y) 3 



;i2) 



where s^ A is the irA invariant mass squared and we take the same functional form (c.f. 
Eq . (|TT1) ) for the 7riVA form factor as for ttNN. The different 7rA charge states are obtained 
from Eq.flTJ) via /^ +A o = (l/2)/ 7r o A + = (l/3)/ ff - A ++ = (1/6) Aa- The vriVA coupling 
constant is defined by: 

(N7r\H int \A) = g„ NA u(p N , s N ) {p% - pi) u a (p A , s A ), (13) 

with u a the Rarita-Schwinger spinor-vector. The value of g^A can be related to the 71NN 
coupling constant via SU(6) symmetry, ^ata = (6V2/5)f-KNN/ m -K ~ H-8 GeV^ 1 (with 
fwNN = ( m 7r/2M) Q-kNn') • Note that in Ref. J|(J (see also f2"2| ) the 7riVA coupling was 
extracted from the width of the decay A — > Ntt, giving a somewhat larger value g^NA — 



15.9 GeV compared with the SU(6) value. However, from numerous studies j|T| of nN 
scattering in the A resonance region, it is known that up to 50% of the A width comes from 
tcN rescattering (for example, through diagrams of the Chew-Low type), so that it would 
be inappropriate to acribe the entire width to the tree level process in determining g n NA- 
We expect the SU(6) value for the 7rAA coupling to be accurate to within ~ 10-20%. 

In Fig. 6 we show the 7rA distribution function as a function of y, compared with the ttN 
distribution (0). The latter is calculated with a form factor cut-off, A n N, of 1 GeV, while 
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the former has A^a = 1-3 GeV, which gives (n) n A = 11%. The 7rA distribution is somewhat 
broader than the irN, peaking at slightly larger y (0.3 c.f. 0.25) and having more strength 
at larger y, so that relatively more of the ttN contribution may be canceled at larger x than 
at smaller x. (Note that in covariant approaches with a t-dependent dipole form factor the 
it A distribution is softer than the ttN, so that there the A plays a negligible role at large 
x.) The cancellation of the d excess with the inclusion of ttA states can be seen explicitly in 
Fig. 7, where the (positive) ttN and (negative) ttA contributions are shown (dashed lines) 
for A„-jv = 1.5 GeV and A^a = 1.3 GeV, together with the sum (solid). In Fig. 7(a) the nA 
brings the difference d — u closer to the large- x data points, while allowing for a reasonable 
fit at smaller x. However, the cancellation is still not sufficient to produce a downturn in 
the ratio at large x, as the E866 data appear to prefer, Fig. 7(b). 

More cancellation can be achieved by either increasing the 7rA contribution, or decreasing 
the ttN contribution. Either is acceptable in the model, as long as the resulting distributions 
do not contradict other observables, such as the total d + u distribution, which should serve 
as an absolute upper limit on the strength of the form factor ||. In Fig. 8(a) we show the 
contributions to the sum x(d+u) from the irN and 7rA components with A^n = 1.5 GeV and 



A^a = 1.3 GeV, compared with the CTEQ4 |~5| and MRS98 [16 parameterizations. While 
at small x the calculated distributions lie safely below the parameterization (the difference 
is made up by the perturbatively generated g — > qq antiquark distributions), at large x 
the pion cloud already saturates the total sea with these cut-offs — although one should 
add a cautionary note that the antiquark distribution at large x is not determined very 
precisely. For softer combinations of form factors, namely A n N = 1 GeV, A^a = 1-3 GeV 
and A^-at = A^a = 1 GeV, the total non-perturbative antiquark sea in Fig. 8(b) is below the 
empirical parameterizations in both cases. 

Therefore the only way to obtain a smaller d excess at large x and still be consistent with 
the total antiquark distribution is to reduce the tcN component, having a cut-off smaller 
than for the tcNA vertex. It was argued in Ref. [I7| that the 7rA^A form factor should be 
softer than the ttNN, based on the observation that the Ml transition form factor was 
softer for 'yNA than for jNN. However, there is no clear connection between these form 
factors, and hence no compelling reason why the nNA form factor cannot be harder than 
that for ttNN. Indeed, a comparison of the axial form factors for the nucleon and for the 
N-A transition strongly favor an N-A axial form factor that is significantly harder than 
that of the nucleon. In fact, the former is best fit by a 1.3 GeV dipole, while the latter by a 
1.02 GeV dipole parameterization [ 4"2"| . Within the framework of PC AC these form factors 



are directly related to the corresponding form factors for pion emission or absorption [43 



In Fig. 9 we show the difference and ratio of the d and u distributions calculated with 
the softer nNN form factor, A n ^ = 1 GeV, and A^a = 1.3 GeV. The excess at large x 
now is largely canceled by the 7rA. However, the smaller nN contribution means that the 
asymmetry is underestimated in the intermediate x range, x ^ 0.2. 

Based on these results it would appear difficult to obtain a quantitative description, 
within the pion cloud model, of both the ratio and the difference of the d and u distributions, 
together with the total sea. One needs to consider, therefore, the possibility that other 
mechanisms may at least be partly responsible for the discrepancy. Additional meson- 
baryon components, such as the pN, could be included in extended versions of the meson 
cloud model P,|2~4"fl. The pN, however, has a harder y distribution than the ttN, which would 
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lead to an enhancement of asymmetry at large x, in contradiction with the data, even though 
the p structure function may be softer than the 7r. On the other hand, the magnitude of the 
pN contribution is known from previous analyses [||,[24[] to be significant only for very hard 
pNN form factors. For a pNN vertex of similar shape to that used for the irNN vertex, 
A P N — A n N, the contribution from the p is unimportant. Furthermore, the pA contribution 
is far too small to be featured in any subsequent cancellation [|5l, p4| , f27| , if the pA form factor 
is comparable to that for nA. 

Another interesting possibility is that the pion sea itself could be asymmetric. In Eq.(^|) 
only the valence structure of the pion was included, though in principle there could even be 
asymmetric contributions to d — u in the proton from asymmetric d n+ and u w+ distributions 
in a pion. One obvious source of a pion sea asymmetry involves the same physics that is 
responsible for the charge radius of the 7r + , namely the dissociation into virtual tt and p 
mesons, 7r + — > 7i + p° or 7r°p + . The effects of a meson cloud of a pion (irp as well as KK* and 
KK ) on deep- inelastic structure functions were previously investigated in Ref. [44|. (Of 



course the A-isobar again cancels some of this asymmetry through the tt~ , but as with the 
valence pion contributions, the sign of the effect remains.) 

Unfortunately, nothing is known empirically about the pion sea, so that the shape and 
normalization of such an asymmetry can at present only be speculative. Rather than con- 
struct a detailed model of the pion sea involving additional free parameters (meson-meson 
vertex functions), at this stage it is more practical to ask how sensitive could the overall 
d/u asymmetry be to a possible asymmetric sea, one that is consistent with all the known 
phenomenological constraints. To address this question we parameterize the non-singlet part 
of the sea distribution in the pion by a simple form, 

d^-ut^Mx^il-xf. (14) 

From low-energy quark models and Regge theory the exponent a is expected to be around 
and -1/2, respectively, while (3 should be between 5-7 from perturbative QCD arguments 
and from our knowledge of the nucleon sea quark distribution. The normalization of this 
component is unknown and given by the parameter N ' . Allowing for up to a factor 2 
uncertainty in the pion sea (which is related to the uncertainty in the knowledge of the 
gluon distribution in the pion), we set J\f ~ 4, and to be definite take a = and (3 = 5. 

The resulting dj u ratio is shown in Fig. 10 for a ratio of pion sea distributions d^ a : u^ a 
= 2:1 (lower solid) and 4:1 (upper solid), respectively, as well as for a symmetric sea (dashed) 
for A^tv = A w a — 1 GeV. Clearly one gets appreciable enhancement in the low- and 
intermediate-x range, bringing the curves to better agreement with the data, even though 
the ratio is somewhat overestimated at low x for the more asymmetric sea scenario. How 
reasonable this choice of parameters is can only be ascertained by acquiring data on the 
pion structure function at values of x smaller than currently available. The phenomenolog- 
ical consequences of an asymmetric pion sea for Drell-Yan vrA^ and other processes will be 



discussed in more detail elsewhere J45] — see also Ref. j|6| for a discussion of measurements 
which would be sensitive to such an asymmetry. 

Going beyond explanations involving meson clouds, one can also investigate the possibil- 
ity that the bare nucleon itself could be asymmetric with respect to u and d. As suggested 
long ago by Field and Feynman J7|, the Pauli exclusion principle can contribute to the 
asymmetry on the basis of the u and d valence quarks being unequally represented in the 



9 



proton, thereby affecting the likelihood with which qq pairs can be created in different flavor 
channels. In fact, earlier analyses of the NMC data |9,10] suggested that the best agreement 
between theory and experiment could be obtained with the combined effects of pions and 
antisymmetrization, and in the next section we explore this possibility further. 



IV. ANTISYMMETRIZATION 

Although not directly attributable to the exclusion principle, the perturbative effects 
of higher-order quark exchange diagrams on the d — u difference were calculated long ago 
by Ross and Sachrajda [18|]. They found that while they had the correct (positive) sign, 
their magnitude was insignificantly small, since they only arose at order a 2 s . The flavor 
asymmetry of the sea associated with the Pauli principle can therefore only be addressed 
within non-perturbative approaches to parton distributions, as concluded in Attempts 
to calculate the valence distribution of the proton in various quark models began in the 



mid-70s [[47|], when it was realized that the relationship to the QCD improved parton model 
is quite natural at a low scale (below 1 GeV 2 ), where most of the momentum of the nucleon 
resides on its valence quarks [pjfl. This observation has been successfully exploited by the 



Dortmund group [J49|, for example, in constructing phenomenological, valence dominated, 



parameterizations in just this region. 

Bag model calculations of nucleon structure functions have provided some interesting 
insights into the non-perturbative parton distributions ||[50|. For any model in which valence 



quarks are confined by a strong scalar field, the vacuum inside and outside the hadron will 
be different. From the point of view of an external probe, such as the virtual photon in deep 
inelastic scattering, the change in the vacuum structure inside the hadron will appear as 
an intrinsic, non-perturbative sea of qq pairs ||. Because of the Pauli exclusion principle, 
the presence of two valence u quarks, as opposed to a single valence d quark, in the proton 
implies an asymmetry in this non-perturbative sea, so that there is a small excess of dd pairs 
over uu pairs. 

For details of the quantitative calculation of this effect, which is model dependent, we 
refer to the original papers ||. It is enough for us that the shape of d — u was found to be 
similar to that of the usual sea quark distributions and the normalization, / dx(d — u), less 
than 0.25. With this in mind, we parameterize the Pauli contribution by: 

(d- w) Pauh = A Pauli (n + 1)(1 - x) n . (15) 

Because the E866 data implies a softer asymmetry than typical global fits of total sea quark 
distributions would give, as Figs. 2 and 3 illustrate, phenomenologically the power n should 
be ^ 10 rather than the 5-7 that has been common for the total q fits [|TR|nj. (Compare 



also with the original Feynman-Field parameterization which had n = 10 and 7 for u and 
d, respectively.) 

The Pauli effect will produce an excess of d over u over the whole range of x, so that 
it cannot lead to any cancellation of the large-x asymmetry. To be consistent with the 
trend of the large-x data, especially for the d/u ratio, one needs therefore to keep the 
nNN contribution softer than that from ttNA. Taking the nN and ttA contributions 
calculated with A^-jv — 1 GeV and A^a = 1-3 GeV as in Fig. 9 above, we show in Fig. 11 the 
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combined effects of pions and the Pauli effect. For the latter the exponent n = 14, and the 
normalization is A Paub pa 7%, which is at the lower end of the expected scale but consistent 
with the bag model calculations M. Together with the integrated asymmetry from pions, 
A 71 " pa 0.05, the combined value A = A 1 + A Pauh pa 0.12 is in quite reasonable agreement 
with the experimental result, 0.100 from E866 and 0.148 from NMC. While the quality of 
the fit in Fig. 11 is quite good, it would be further improved (see Fig. 12) if one were to use 
the softer pion structure function, q* ~ (1 — x) 2 , as suggested by perturbative QCD 



Before leaving this discussion of antisymmetrization we should also mention the calcu- 
lations of Donoghue and Golowich J51[| , and more recently Steffens and Thomas [Q, of the 
one-gluon-exchange corrections to the 3-quark proton wave function. Considering all possible 
permutations of the 5-quark wave function allowed by Fermi statistics, the antisymmetriza- 
tion of the qq pair split off from the emitted gluon with the quarks in the nucleon ground 
state, it turns out that there are in fact more diagrams for u quarks than d quarks. This pe- 
culiarity results in there actually being an excess of u quarks over d, albeit a very small one. 
On the other hand, one should note that this calculation considered just the perturbative 
contribution, while the Signal-Thomas effect || is a totally non-perturbative phenomenon, 
including all possible non-perturbative interactions between the produced quark (or anti- 
quark) and the confining mean field of the proton. Steffens and Thomas also investigated 
the effects of antisymmetrization between qq pairs arising from one-pion loops with the three 



quarks in the nucleon ground state |p2| , although here again the effects were found to be 
quite small compared with the antisymmetrization for the bare nucleon state, and from the 
pion cloud contribution discussed in Section III. 



V. CONCLUSION 

We have, for the first time, at our disposal important new data which map out the 
x-dependence of the asymmetry of the light antiquark sea. Most importantly, the E866 
Drell-Yan results confirm the earlier observations that the d and u content of the proton 
is not symmetric. One of the more interesting new features of the data is the relatively 
fast downturn in the d/u ratio beyond x ~ 0.15, which drops rapidly back to unity by 
x ~ 0.3. Taken at face value, this would appear to provide a challenge to models in which 
the asymmetry is assumed to arise solely from the pion cloud of the nucleon, and in turn 
leads us to consider a richer and more complex structure of the non-perturbative sea in 
which several mechanisms may give competing contributions. 

The evidence from the large- x data indicates that a ttA component in the nucleon wave 
function may be necessary, one which is harder in momentum space than the irN component. 
Such a distribution arises naturally in the infinite momentum frame formulation of the pion 
cloud, unlike in earlier covariant approaches using t-dependent form factors where it was 
softer than the nN component and hence played no role at large x. Consistency with 
data for the sum of d and u at x > 0.2 requires that both the ttNN and ttNA form 
factors be relatively soft, making it difficult to avoid underestimating the E866 asymmetry 
at intermediate x, and leaving room for other effects, such as the Pauli exclusion principle, to 
make up the difference. Along the lines of previous estimates of the Pauli effect, we find the 
contribution to the d — u difference from antisymmetrization to be significant in magnitude, 
and particularly important at small x. Our final results suggest that the best description 
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of the E866 data is indeed that in which pions and antisymmetrization play roughly equal 
roles — consistent with the findings of the earlier analysis [|l(J of the NMC data for F^ — F^. 

In conclusion, we note that it would be helpful to have more data at large x, where the 
error bars are largest, to verify the downward trend of d — u, and to further explore the 
possible discrepancy between the Fermilab and CERN data. 
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FIG. 1. Ratio of the antiquark distribution in a nucleon bound in the deuteron to that in the 
free nucleon, for q ~ (1 — x) n , with n = 5, 7 and 10. 
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FIG. 2. The x dependence of the d/u ratio from the E866 p|] (filled circles) and NA51 M (open 
circle) experiments, compared with the CTEQ4 |l5| and MRS98 parameterizations. Note that 
the MRS98 parameterization included the E866 data in their fits, while CTEQ4 predates the 
experiment. 
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FIG. 3. Comparison of d — u from the E866 experiment [I]] with the values extracted from 
the NMC measurement of the proton-neutron structure function difference Q, using the CTEQ4 



[15] and MRS98 [16| parameterizations for the valence quark distributions. Also shown are the 
parameterizations of d — u from CTEQ4 (dashed) and MRS98 (dotted). 
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FIG. 7. Contributions from the irN and irA components (dashed) and the combined effect 
(solid) to the (a) d — u difference and (b) d/u ratio. The cut-off masses are A^n = 1.5 GeV and 
A^A = 1.3 GeV. 
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FIG. 8. Total x(d + u) distribution (a) from the ttN and irA components (dashed), with 
Ajttv = 1-5 GeV, A„-a = 1.3 GeV, and the total (solid), (b) the total contribution for 
A-n-jv = 1.5 GeV, A^a = 1-3 GeV (largest curve), A nN = 1 GeV, A^a = 1-3 GeV (middle), 
and A^iv = A„a = 1 GeV (smallest). The theoretical curves are compared with the CTEQ4 Hf 
and MRS98 [16| global parameterizations (dotted). 
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FIG. 10. Effect of an asymmetric pion sea on the d/u ratio. The dashed curve represents the 
ratio for a symmetric pion sea with A w n = A„-a = 1 GeV, while the solid curves have asymmetric 



seas in the ratio d^ ea 



u" 



2 : 1 (lower curve) and 4 : 1 (upper curve). 
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FIG. 11. Contributions from pions with A^tv = 1 GeV and A^a = 1-3 GeV (dashed) and from 
antisymmetrization (dotted) to the (a) d — u difference and (b) d/u ratio, and the combined effect 
(solid). 



22 



2 



I 

\ 
1^ 



0.5 





i 


i 










> 






i 

■ 





0.1 0.2 0.3 

JO 

FIG. 12. As in Fig. 11(b), but with an extra power of (1 — x) in the pion structure function, 
according to Ref . pl| . 



23 



